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Abundance determination from global emission-line SDSS spectra: 
exploring objects with high N/O ratios 
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ABSTRACT 

We have compared the oxygen and nitrogen abundances derived from global emission-line 
SDSS spectra of galaxies using (1) the T e method and (2) two recent strong line calibrations: 
the ON and NS calibrations. Using the T e method, anomously high N/O abundances ratios 
have been found in some SDSS galaxies. To investigate this, we have Monte Carlo simulated 
the global spectra of composite nebulae by a mix of spectra of individual components, based 
on spectra of well-studied H n regions in nearby galaxies. We found that the T e method results 
in an underestimated oxygen abundance (and hence in an overestimated nitrogen-to-oxygen 
ratio) if H n regions with different physical properties contribute to the global spectrum of 
composite nebulae. This effect is somewhat similar to the small-scale temperature fluctua- 
tions in Hn regions discussed by Peimbert. Our work thus suggests that the high r e -based 
N/O abundances ratios found in SDSS galaxies may not be real. However, such an effect is 
not expected to be present in dwarf galaxies since they have generally an uniform chemical 
composition. The ON and NS calibrations give O and N abundances in composite nebulae 
which agree with the mean luminosity-weighted abundances of their components to within 
-0.2 dex. 

Key words: galaxies: abundances - ISM: abundances - H n regions 



,H ,1 INTRODUCTION 



Metallicities play a key role in many studies of galaxies. Gas- 
phase oxygen and nitrogen abundances are broadly used to esti- 
mate these metallicities. It is believed (e.g. IStasi hska 200(f) that 
emission lines in spectra of Hn regions are the most powerful 
indicators of the chemical composition of galaxies, both in the 
low- and intermediate-redshift universe. The spectra of a large 
number of individual Hn regions in nearby spiral and irregu- 
lar ga l axies have now been obtained ( McCall, Rvbski & Shields 
19851: IZaritskv, Kennicutt & Huchral 1 1994 Ivan Zee et al. I 1 199% 
van Zee & Havn es I 20061; iBresolin. Kennicutt. & Garnett I 1 19991 ; 
Bresolin et a l. 2005, 2009, among many others). These spectro- 



due to the completion of several large spectral surveys such as 
the Sloan Digital Sky Survey (SDSS, lYork et al. I l200d) . Mea- 
surements of emission lines in SDSS spectra have been used for 
abundance determinations in a number of studies JKniazev et al. I 
20041: llzotovetal. Il200fj| : iTremonti et aTll2004 lAsari et al. 1120071 : 



Thuan et al. 120101 among others). The auroral lines are measurable 



in a relativel y large number of SDSS galaxies jKniazev et al. 120041 : 
2006 ) that provide the possibility to obtain 7Vbased 



in a relatively 
llzotov et alT : 



scopic measurements provide basis for investigations of metal- 
licity properties (such as radia l abundance gradients, cen- 
tral metallicities, etc) of galaxies ([Za ritskv. Kennicutt & Huchra 
19941: Ivan Zee et al. I \l99$: IBresolin, Kennicutt, & Garnett 1 1 199*91 ; 
Pilyugin. Vflchez. & Contini Il2004bl . among others). 



Kennicutt ] h992l) pioneered another method of spectral in- 
vestigation of galaxies, considering the global, i.e. spatially un- 
resolved, spectrophotometry of a sample of 90 galaxies spanning 
the entire Hubble sequence. In recent years, the number of avail- 
able spectra of emission-line nebulae has increased dramatically 



abundances for SDSS galaxies. 

The SDSS spectra are obtained through 3-arcsec diameter 
fibers. At a redshift of z = 0.025, the projected aperture diame- 
ter is ~ 1.5 kpc, while it is ~ 10 kpc at a redshift of z = 0.17. This 
suggests that SDSS spectra of distant galaxies are closer to global 
spectra of galaxies, i.e. they are the integrated spectra of multiple 
rather than just one H n region. Therefore the meaning of 7Vbased 
abundances in SDSS galaxies is unclear. If a single giant Hn re- 
gion, excited by a single star cluster, is responsible for the global 
SDSS spectrum then one can expect that the 7Vbased abundance 
to be a good estimate of the true one. Such a situation occurs only 
in nearby SDSS galaxies and in some distant SDSS galaxies where 
a single supergiant Hn region, resulting from a strong starburst, 
makes a dominant contribution to the global spectrum. One can ex- 
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pect, however, that in a majority of distant SDSS galaxies, multiple 
individual H n region s contribute t o the global spectrum. 

lErcolano et ail {2007, 2010) have investigated the effect of 
multiple ionization sources in Hn regions on the total elemental 
abundances derived from the analysis of collisionally excited emis- 
sion lines. They considered the case of an ionizing set composed 
of two stellar populations with masses of 37 Mq and 56 Mq. They 
have shown that the temperature structures of models with centrally 
concentrated ionizing sources can be quite different from those of 
models where the ionizing sources are randomly distributed within 
the volume, with generally non-overlapping Stromgren spheres. 
Since the SDSS composite nebula can be ionised by sources of 
different temperatures then the meaning of the measured electron 
temperature in distant SDSS galaxies and consequently, that of the 
7Vbased abundance is not evident. 

This problem is in some sense similar to the one concerning 
the validity of the 7Vbased abu ndances in H n r egions with small- 
scale temperature fluctuations I Peimb ertll 19671) or/and with large 
variations o f the electron temperature across the nebula JStasinskal 
Il978l l2005h - According to Peimbert, since the line fluxes in the 
spectrum of such a nebula are weighted in favor of the hot regions, 
the electron temperature estimated from the auroral-to-nebular 
lines ratio will be overestimated and the abundance underestimated. 
While the existence of such temperature fluctuations in individual 
H ii regions is still debated, temperature variations of individual H n 
regions within a composite nebula are to be perfectly natural. Thus, 
such an effect is probably at work in the distant SDSS galaxies. 

When the electron temperature of an extragalactic Hn re- 
gion cannot be measured, then its location in some emission-line 
diagrams is used for estimating its oxygen abundance. This ap- 
proach to a b und ance deter mination i n Hn regions, proposed by 
IPagel et al. I i 19791) and lAlloin et all dl979l) . is usually referred 
to as the "strong line method". Numerous relations have been 
proposed to convert metallicity-sensitive emission-line ratios into 
metallicity or temperature est i mates ( e.g. iDopita & Evans I 
Zaritskv. Kennicutt & Huchra I ["l994l; IVilchez & Esteban 
Pilvuginl |2000| l200ll; iDenicolo, Terlevich & Terlevich I 120021 
Pettini & Pageil |2004|; iTremonti et all |2004 IPilvugin & Thuanl 
2001 iLianget alJ 12004 IStasinskal 120061; jBresolinl l2007h 
Perez-Montero & Contini I 120091 : iThuan et al. I l20ld) . Although 



1986 



1996 



SDSS spectra of distant galaxies are closer to global galaxy spec- 
tra than to spectra of individual H 11 regions, abundances in distant 
SDSS galaxies have been estimated using the strong line meth- 
ods de veloped for abundance determination in individual Hn re- 
gions dTremonti et al. I [20041: lErb et all 120061 ; lAsari et al. I [20071 ; 



Thuan et al. 



2010, among others). This, despite the warning of 



Stasinska ~ j2010h that the strong line methods should be used only 



for nebulae having the same structural properties as those of the 
calibration samples. 

However, there is a reason to expect the calibrations may pro- 
vide more robust abundances in composite nebulae than the T e 
method. For definiteness, we will discuss the case of oxygen abun- 
dances. The auroral and nebular lines originate in transitions from 
levels that differ significantly in energy, with the levels that give rise 
to the auroral line being at higher energies than those for nebular 
lines. Therefore, the emissivity of the auroral line depends much 
more strongly on the electron temperature in the nebula than that 
of the nebular line. Hence, the auroral and nebular line fluxes in 
the integrated spectra of composite nebulae are weighted in differ- 
ent ways: the hot regions will give relatively more weight to the 
auroral line than to the nebular lines. This means the nebular and 
auroral lines in the spectra of composite nebulae correspond to dif- 



ferent temperatures, i.e. they are not self-consistent. The electron 
temperature estimated from the auroral-to-nebular lines ratio will 
then be overestimated and the abundanc e will be underestimated 
jPeimbert|[l967l : IStasinskalll978l l2005h . In the case of the em- 
pirical calibrations considered here (the ON and NS calibrations), 
combinations of strong lines for different ions serve as abundance 
and electron temperature indicators. The differences between the 
energies of parent levels for these nebular lines are smaller than 
those between the energies of parent levels for the auroral and neb- 
ular lines. Then one can expect that, in the strong line method, the 
different regions will give roughly similar relative contributions to 
the nebular line fluxes, i.e. they will all correspond to more or less 
similar electron temperatures. This, in turn, will result in more ro- 
bust global abundances estimates. 

The goal of the present study is to examine to what extent 
abundances for distant SDSS galaxies derived in different ways (by 
the classic T e method and through recent calibrations) agree or dis- 
agree. The abundances of a sample of SDSS galaxies are examined 
in Section 2. In Section 3, we compare the data with artificial global 
spectra of composite nebulae generated by Monte Carlo techniques 
as mixes of spectra of individual components, using spectra of real 
Hn regions in nearby galaxies with measured electron tempera- 
tures. The meaning of abundances derived in composite nebulae 
in different ways and their locations in the O/H - N/O diagram are 
examined. A discussion of the results is given in Section 4. Section 
5 presents the conclusions. 

Throughout the paper, we will be using the following standard 
notations for the line intensities: 

Rl = I{OU]Aiin+AH29l 'At/J, 
Nl = /[«//] / 16548+/16584/-/h/J, 



<[OI1I],14959+,(5007 



The electron temperatures will be given in units of 10 K. 



2 OXYGEN AND NITROGEN ABUNDANCES 
2.1 Galaxy samples 

2.1.1 SDSS galaxies 

In our previous work JPilvugin et al~ll2010ah . we have constructed 
a sample of galaxies by selecting from Data Release 6 of the SDSS 
the galaxy spectra which satisfy the following criteria: (1) both 
[Om]^4363 and [O n]/17320+/17330 auroral lines are detected; (2) 
the spectra have smooth line profiles. Particular attention was paid 
to the two auroral lines mentioned above since the accuracy of the 
electron temperature determination depends mainly on those emis- 
sion lines; (3) the emission lines do not have a broad component. 
The line intensities in these spectra have been measured by fitting 
every line with a Gaussian profile a n d de-reddened in the w ay de- 
scribed in lPilvugin&Thuanl j2007l) : |Pirvugin et al. I feoiOal) . 

The wavelength range of the SDSS spectra is 3800 - 
9300A. Hence, for nearby galaxies with redshift z ;£ 0.02, the 
[On]/13727+/13729 emission line is outside that range. The ab- 
sence of this line prevents determination of the oxygen abundance 
through the standard version of the T e method. We thus exclude 
the most nearby galaxies. All galaxies in our sample have redshifts 
larger than ~ 0.023, i.e. they are more distant than -100 Mpc. We 
thus obtain a total sample of 281 SDSS galaxies. The upper panel 
of Fig. Q] shows the positions of the SDSS galaxies of our sample 
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Figure 1. The logif3 - logi?2 diagram. The filled dark (black in the color 
version) circles show the objects from our SDSS sample (upper panel), the 
sample of calibration Hn regions (middle panel), and the Green Pea galax- 
ies (lower panel). The filled grey (light-blue in the color version) circles in 
each panel are objects from the sample of Hn regions in nearby galaxies. 
(A color version of this figure is available in the online version.) 

in the standard if 3 - R 2 diagram by filled dark (black in the color 
version) circles. 



individual Hn regions in nearby spiral and irregular galaxies. Their 
sample consists of 1121 data points and will be used to outline 
the area occupied by Hn regions of nearby galaxies in the if 3 - 
if 2 diagram. These are shown by filled grey (light-blue in the color 
version) circles in the if 3 - if 2 diagram. Fig.Q]shows that the SDSS 
galaxies are located in a relatively small part of the area occupied 
by H 11 regions in nearby galaxies. 



2. 1.3 Calibration H n regions 

The [Om]/14363 or/and [Nn]/15755 auroral lines are detected in 
many H 11 regions i n nearby galaxies. A com pilation of such H 11 re- 
gions was made bv lPilvugin et al. IfeOlObl) and used as calibration 
data points to derive the relations which give the oxygen and ni- 
trogen abundances and electron temperature in terms of the strong 
emiss ion line fluxes. After the exclusion of a few "peculiar" ob- 
jects JPilvugin et aLll2010bl) . we end up with 112 data points. The 
spectra of this sample of H 11 regions in nearby spiral and irregu- 
lar galaxies with detected auroral lines will be used as templates to 
simulate artificial spectra of composite nebulae. The middle panel 
of Fig. [T] shows the positions of the calibration Hn regions in the 
if 3 - R2 diagram as filled dark (black) circles. The filled grey (light- 
blue) circles are the same as in the upper panel. 



2.1.4 Green Pea galaxies 

The G reen Pea galaxies have be en extracted from the SDSS sam- 
ple by ICardamone et al. I 12009). These galaxies have a compact 
appeara nce and are charac t erised by distinct green colours in gri 
images. ICardamone et al. ] J2009h found that the Green Peas are 
low-mass galaxies (M ~ 10 85 - 10'°Mq) with high star for- 
mat ion rates (~ I OMq yr~') with metallicities 12 +log(0/H) ~ 
8.7. Izotov et al. ] d201lh have argued that they form a subset of 
a larger population of luminous compac t star-forming galaxies. 
Amonn, Perez-Montero & Vflchez (2010) have derived r e -based 
oxygen abundances in these galaxies and have found that the Green 
Pea galaxies show metallicities 7.7 < 12 + log(0/H) < 8.4, with a 
mean value of 8.05 ± 0.14. They have al so found that some Green 
Pea galaxies display enhanced N/O ratios. llzotov et al. 1) have 
found that the oxygen abundances 12 + log(0/H) in luminous com- 
pact star-forming galaxies are in the range 7.6 - 8.4. They found 
no appreciable difference in element abundances between luminous 
compact star-forming galaxies and local blue compact dwarf galax- 
ies, implying a similar chemical enrichment history. 

From the sample of Green Pea galaxies of Carda mone et al. I 

d2009h we have selected those with detectable [O m]/14363 auro- 
ral lines. The line fluxes in the SDSS spectra were measured with 
irafQ Th e measured line fluxes we r e dereddened in the wa y de- 
scribed in Pilyugin & Thuan ] j2007h ; |Pilvuginetal.l d2010ij) . The 
lower panel of Fig. [T] shows the positions of the Green Pea galax- 
ies in the log if 3 - log if 2 diagram (filled dark (black) circles). The 
filled grey (light-blue) circles are the same as in the upper panel. 
Note that the Green Pea galaxies overlap with the area occupied by 
our SDSS sample in the if 3 - if 2 diagram. 



2.1.2 Hn regions in nearby galaxies 

IPilvugin. Vflchez, & Contini I d2004bl) have made a compilation of 
a large amount of strong emission line measurements in spectra of 



iraf is distributed by National Optical Astronomical Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



© 0000 RAS, MNRAS 000, 000-000 



4 L.S.Pilyugin et al. 



2.2 Abundance determination 

Line fluxes are converted to electron temperatures and ion abun- 
dances within using the standard Hn region model with two dis- 
tinct temperature zones within the nebula: the electron tempera- 
tures t 2 within the + zone and within the ++ zone. Electron 
temperatures and ion abundances fro m line fluxes were derived as 
according to Pilvugi n et al. | l l2010bl) . The ratio of the nebular to 
auroral oxygen line intensities [Om]/i4959 + /15007/[O m]/l4363 is 
used for the r 3 determination and the ratio of the nebular to auroral 
nitrogen line intensities [N ii]/16548 + /16584/[N n]/i5755 is used for 
the ti determination. 

When only one of the electron temperatures is known, it is 
common practice to estimate the other temperature by using a rela- 
tion between t 2 and ts . We have adopted the following t2~t3 relation 

t 2 = 0.672/ 3 +0.314. (1) 

This relation derived by Pilvugin et al. ( 2009) i s very similar to 
the wi dely used one pro posed by [Campbell. Terlevich & Melnickl 
dl986h and confirmed bv lGarnettl j 19921) . In general, using the clas- 
sic T e method, one can derive the oxygen abundance in each Hn 
region in several ways: (0/H), 3 abundances can be determined with 
the measured tj, t 2 being then estimated from the t2-t3 relation; or 
(0/H), 2 abundances can be found with a measured t 2 , h being then 
estimated from the h-t^ relation. If the line-flux measurements are 
accurate enough, the oxygen abundances derived in these two ways 
should agree. The (O/H),, abundance will be considered here and 
will be referred to as (0/H) Te hereafter. 

We also estimate the oxygen and nitrogen abundances us- 
ing two recent empirical calibrations. The ON-calibration rela- 
tions give the oxygen and nitrogen abundances and electron tem- 
perature in terms of the fluxes of the strong emission lines ++ , 
+ , and N + . It has been derived using spectra of Hn regions 
with well-measured ele ctron temperatures as calibration datapoints 
jPilvugin et al. I [201 Obi) . The oxygen and nitrogen abundances es- 
timated using the ON calibration will be referred to below as 
(0/H)on and (N/H) 0JV . 

In addition to the direct method and the ON-calibration, the 
NS-calibration is used as well to estimate oxygen and nitrogen 
abundances. These relations give abundances and electron tempera- 
tures in terms of the fluxes in the strong emission lines of ++ , N + , 
and S + , and have also been derived using spectra of Hn regions 
with well-measured electron temperatures as calibration datapoints 
jPilvugin & Mattsson 1201 lh . The oxygen and nitrogen abundances 
estimated using the NS calibration will be referred to below as 
(O/HXvs and (N/H) ws . 



2.3 The O/H - N/O diagram 

Fig- IU shows a standard O/H - N/O diagram for the calibration H n 
regions. Using these data, we have obtained the following linear 
relations between log N/O and log O/H 

log (N/O) = -1.493 

for 12 + log (O/H) < 8.14, 

= 1.489 x (12 + log (O/H))- 13.613 
for 12 + log (O/H) > 8.14. 



(2) 



Several H n regions show large (>0.25 dex) abundance deviations. 
These have been excluded in deriving the final relation. In Fig. [2] 
the data points used in the derivation of the final relation (102 out 
of 112) are shown by filled circles, while the rejected data points 
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Figure 2. The O/H - N/O diagram for the calibration H n regions. The solid 
(black) line shows a linear fit to the N/O - O/H relation (Eq.FJJ. The filled 
circles are data points used in deriving this relation, the open circles are 
rejected data points. The dashed (red) line shows an alternative quadratic fit 
to the N/O - O/H relation. 



are shown by open circles. The solid (black) line shows the derived 
relation. An alternate quadratic relation between log N/O and log 
O/H has also been derived: 



log (N/O) = - 1.39 + 0.39 (Z-Z ) + 0.30 (Z-Z ) 2 
for Z < 8.23, 

= -1.39 + 1.24 (Z-Z ) + 1.63 (Z-Zo) 2 
for Z > 8.23, 



(3) 



where Z = 12 + log(0/H) and Z = 8.23 (see the dashed (red) line in 
Fig. |2j. The two relations are similar for metallicities 12 + log(0/H) 
< 8.5. We will here only consider objects in this metallicity range. 
Using the N/O - O/H relationship, one can estimate the N/O ratio 
which corresponds to a given oxygen abundance. The nitrogen-to- 
oxygen abundance ratio obtained in this way will be referred to 
as (N/0) R£L below. For definiteness, we will use the N/O - O/H 
relationship given by Eq. [2] 

We have derived (O/H)^, (0/H) 0iv , (O/H)^ oxygen abun- 
dances and (N/H)^,, (N/H)o,v, and (N/H)^ nitrogen abundances 
for all calibration H n regions, SDSS and Green Pea galaxies. Fig.|3] 
shows the comparison between O/H - N/O diagrams for abun- 
dances computed in different ways. 

The upper panel of Fig. [3] shows the O/H - N/O diagram 
for SDSS objects (filled dark (black) circles), Green Pea galax- 
ies (dark (red) plus signs), and calibration H n regions (filled grey 
(light-blue) circles) for abundances obtained with the T e method. 
The upper panel shows that some SDSS objects deviate signif- 
icantly from the general trend defined by the calibration Hn 
regions: they are shifted towards higher N/O abundances ratio 
or/and towards lower O/H abundances. The same shift is observed 
for some Green Pea galaxies, in agreeme nt with the findings of 
lAmorfn. Perez-Montero & Vflchez I d2010)) . 

The middle panel of Fig. [3] shows the same O/H - N/O dia- 
gram, but for the case where abundances are obtained with the ON 
calibration, and shows that the SDSS and Green Pea galaxies follow 
the trend defined by the calibration H n regions, without large de- 
viations. The lower panel of Fig.[3]shows theNS calibration-based 
abundances. Again, the SDSS and Green Pea galaxies follow the 
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Figure 3. The 0/H - N/O diagrams for SDSS objects (filled dark (black) 
circles), for the Green Pea galaxies (dark (red) plus signs), and for calibra- 
tion H ii regions (filled grey (light-blue) circles) with abundances obtained 
with the T e method (upper panel), using the ON calibration (middle panel), 
and the NS calibration (lower panel). (A color version of this figure is avail- 
able in the online version.) 
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Figure 4. The deviation of the N/O ratios in SDSS objects relative the N/O 
- O/H relation as a function of (1) the difference between T e -based and 
ON-calibration-based oxygen abundances (upper panel), (2) of the flux in 
the auroral line [O m]/14363 normalised to the H/J line flux (middle panel), 
and (3) of the equivalent width EW([0 iii]/14363) of the auroral line (lower 
panel). 



trend defined by the calibration H n regions and do not show large 
deviations. 

The upper panel of Fig. [4] shows the deviation of the (N/0)r e 
ratio in the SDSS objects from the N/O - O/H relation as a 
function of the difference between r f -based and ON-calibration- 
based oxygen abundances. There is a clear anti-correlation between 
0og(N/O) r , -log(N/0) fi£i ) and (log(0/H) rf -log(0/H) OAr ) values. 
This anti-correlation suggests the deviation of SDSS objects from 
the general N/O - O/H trend in the upper panel of Fig. [3] is mainly 
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caused by shifts towards lower O/H abundances, i.e. the (0/H) rf 
abundances in these objects appear to be underestimated. 

One may expect the high (N/0)r f ratios in SDSS objects to be 
caused by uncertainties in the flux of the weak auroral [O in]/14363 
line in these objects. To check this hypothesis, we show in the mid- 
dle panel of Fig. [4] the deviation of (N/0)r e ratios (log(N/0)r e - 
log(N/0)/;£L) versus the flux of the [Om]^4363 auroral line, nor- 
malised to the Hy6 flux. The lower panel of Fig.|4]shows the devia- 
tion of the same (N/0)r e ratios as a function of the equivalent width 
EW([Om]^4363) of the [Om]/l4363 auroral line. Examination of 
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Figure 5. The O/H - N/O diagram for artificial composite nebulae and com- 
ponents. Upper panel. Filled dark (black) circles show the r e -based abun- 
dances of components, while the grey (light-blue) filled circles are abun- 
dances of artificial composite nebulae, calculated as the mean luminosity- 
weighted abundance of all components, for case A. Lower panel. The same 
as in the upper panel, but for case B. (A color version of this figure is avail- 
able in the online version.) 



these two panels suggests that the uncertainties in the Om]/i4363 
line flux can indeed be responsible for a significant fraction of high 
(N/0)r e ratios in the SDSS galaxies. However, high (N/0)r e ratios 
are found not only in SDSS objects with a weak [Om]/i4363 line, 
but also in some objects where that line is relatively strong, with 
EW([Om]i4363) > 5 and/or [O m]^4363/H/3 > 0.05. This implies 
that the auroral line flux uncertainties cannot be the only reason for 
the high (N/0)r f ratios found in some SDSS objects. 

We propose here another explanation for the derived high 
(N/O)^ ratios. SDSS galaxies are composite nebulae which con- 
tain a number of H n regions with different physical properties, all 
contributing to the global spectrum. In that case, the meaning of a 
single measured electron temperature for the whole galaxy is not 
well defined and consequently, 7Vbased abundances are ambigu- 
ous. We expect the electron temperature estimated from the auroral- 
to-nebular lines ratio to be overestimated and the derived oxygen 
abundance to be underestimated. In other words, such nebulae will 
ex hibit an eff e ct sim ilar to the temperature fluctuations discussed 
by IPeimbertl d 19671) . However, these temperature variations will 
not occur within a single H n region, but on the large scale, in on- 
going from one H n region to other. These temperature variations 
in an object will cause it to deviate from the general N/O - O/H 
trend. To explore this scenario quantitatively, we use a Monte Carlo 
simulation to construct spectra of composite nebulae (described in 
the next section), and examine how derived abundances in com- 
posite nebulae change, depending on the abundance determination 
method. 



3 GLOBAL ABUNDANCES IN ARTIFICIAL 
COMPOSITE NEBULAE 

3.1 Monte Carlo simulation of global spectra 

To clarify the reason for the discrepancy between the 7Vbased and 
calibration-based abundances in our sample of distant SDSS galax- 
ies, we have simulated their spectra. We have remarked above that 
several different Hn regions can contribute to the global spectra 
of distant SDSS galaxies. We have modelled their globa l spectr a 
by using the approach of IPilvugin. Contini. & Vflchez I {2004a ); 
IPilvugin et ail < l2010ij) . The artificial spectra of SDSS galaxies 
have been computed as a mix of spectra of individual components, 
which we take to be spectra of real Hn regions in nearby galax- 
ies with measured electron temperatures, taken from the calibration 
H ii region sample. 

The auroral line [Om]/14363 is usually detected in spectra of 
low-metallicity calibration Hn regions and t 3 can then be deter- 
mined, while the auroral line [Nn]/15755 is detected in spectra of 
high-metallicity calibration Hn regions, resulting in the determina- 
tion of ?2- When [Nii]/15755 is detected in the spectrum of an Hn 
region, then the intensity of [O m]/14363 is estimated in the follow- 
ing manner. We calculate t 3 from the measured t 2 using the t 2 - t 3 
relation. We then estimate the line flux of [O m]/14363 correspond- 
ing to the obtained value of . 

Using the calibration H n regions as components, we have per- 
formed a variety of Monte Carlo simulations, producing in each run 
an artificial spectrum of a composite nebula. The calibration Hn 
regions are not distributed uniformly within the considered metal- 
licity range. We have thus divided the metallicity range in bins of 
A(log(0/H)) = 0.05 and selected one calibration H n region within 
each bin according to the following two criteria: 

(i) The average value of the difference between the calibration- 
based and T e - based abundances is minimum. Both oxygen and 
nitrogen abundances and both ON and NS calibrations were taken 
into account in determining the average value of the difference. 

(ii) The deviation of N/O abundances ratio from the N/O - O/H 
relation is not in excess of 0.15 dex. 

Fig.[5]shows the positions of the selected calibration H n regions in 
the O/H - N/O diagram by dark (black) filled circles. 

At a fixed 12+log(O/H) , we have simulated 100 global spec- 
tra, using the spectra of selected calibration H n regions with oxy- 
gen abundances within 12+log(0/H)o ± A(log(0/H)). The contri- 
bution from each component to a given global spectrum of artificial 
composite H n regions is defined by the value w s , where wj are ran- 
dom numbers between -1 and 1. To take into account that some 
components may not contribute to a given spectrum we choose the 
interval of random numbers between -1 and 1 instead of the stan- 
dard interval between and 1 and when wj < we set it to 0, i.e. 
the component does not make contribution to the global spectrum. 
The component (calibration H n region) line fluxes are given on a 
scale where F(Hf}) = 1, therefore we can consider that the value w- s 
defines H/? flux (or luminosity) of j component in the given variant 
of the global spectra. Then the total F(H/3) flux of the composite 
nebula is given by the expression: 



(4) 



where n is the number of components in the composite nebula. 
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Figure 6. The O/H - N/O diagrams for the artificial composite nebulae for case A (left panels) and for case B (right panels). The grey (light-blue) filled circles 
show the abundances derived in different ways from artificial global spectra. The dark (black) filled circles are 7" e -based abundances in the components. (A 
color version of this figure is available in the online version.) 



The flux F(X^ k ) of the composite nebula in the line X^ k is 
given by: 

= >lL — • (5) 
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We have computed the intensity of the H/3, [On]/13727, 
[Om]/l4363, [Om]A4959, [OmJ/15007, [Nn]A6548, [Nn]^6584, 
[Sn]/i6717, and [Sn]/i6731 lines for the each artificial spectrum 
of composite nebulae. 

The mean oxygen abundance of a composite Hn region, 
weighted by the H/3 line luminosity, can be determined as 
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Z Wj(0/H)j 

(0/H),„ cfl „ = J -^— . (6) 

E Wj 

Similarly, the mean nitrogen abundance in a composite H n region 
is given by 

2>j(N/H)j 

(N/H),„ cfl „ = J -^— . (7) 

2>j 

For each computed spectrum, we have estimated t 3 , (0/H) mfn „, 
(N/H),„ ra ,„ (0/H) Te> (N/H) r „ (O/HW, (N/H) ws , (0/H) 0(V and 
(N/H) 0lV abundances. This procedure was repeated for a variety of 
12+log(O/H) values . 

3.2 Global abundances in artificial composite nebulae 

Here, we discuss two cases, represented by two sets of Monte Carlo 
simulations: 

A. We simulated 100 global spectra for each 12+log(O/H) , 
using components with oxygen abundances within 12+log(0/H)o 
± 0.15, and varying 12+log(O/H) from 7.5 to 8.5 in steps of 0.05 
dex. The grey (light-blue) filled circles in the upper panel of Fig. [5] 
show the mean abundances of the artificial composite nebulae for 
case A. The abundances of each component has been weighted by 
the corresponding H/3-luminosities. 

B. We simulated 100 global spectra for each 12+log(0/H)o, using 
components with oxygen abundances within 12+log(O/H) ± 0.35. 
The grey (light-blue) filled circles in the lower panel of Fig.[5]show 
the mean abundances in artificial composite nebulae for case B. 

The N/O - O/H diagrams for case A with abundances ob- 
tained in different ways from the global spectra are shown in the 
left panels of Fig. [6] by grey (light-blue) filled circles. The up- 
per panel shows the 7Vbased abundances, the middle panel the 
ON-calibration-based abundances, and the lower panel the NS- 
calibration-based abundances in the artificial composite nebulae. 
The dark (black) filled circles in each panel are the 7Vbased abun- 
dances of each component. The upper left panel shows that a sig- 
nificant fraction of the artificial composite nebulae with T e - based 
abundances are shifted from the O/H - N/O relation towards higher 
N/O ratios or/and towards lower oxygen abundances. 

The left upper panel of Fig. [7j shows the difference between 
nitrogen-to-oxygen ratios (N/0) rf and (N/0) m( , a „ in the artificial 
composite nebulae as a function of the difference between oxygen 
abundances (0/H) r< , and (0/H)„ Kn „ for case A by the grey (light- 
blue) filled circles. The dark (black) filled circles show the com- 
ponents. The right upper panel of Fig. [7j shows the deviation of 
(N/O)^ ratios in the artificial composite nebulae from the N/O - 
O/H relation (the linear fit discussed previously) as a function of 
the difference between the oxygen abundance derived by the T e 
method from global spectra and that obtained as the H/Muminosity- 
weighted mean abundances of components for case A. 

The N/O - O/H diagram of 7Vbased abundances for artifi- 
cial composite nebulae for case B is shown in the upper right panel 
of Fig. [6] by the dark (black) filled circles. The grey (light-blue) 
filled circles are 7Vbased abundances of the components. Again, a 
fraction of artificial composite nebulae with 7Vbased abundances 



are shifted from the O/H - N/O relation towards higher N/O ra- 
tios or/and towards lower oxygen abundances. The grey (light-blue) 
filled circles in the lower left panel of Fig. |7J shows the difference 
between nitrogen-to-oxygen ratios (N/0)r t and (N/0) mfM „ in arti- 
ficial composite nebulae as a function of the difference between 
oxygen abundances (O/H)^ and (OfH) mea „ for case B. The dark 
(black) filled circles show the individual components. The lower 
right panel of Fig. UJ shows the deviations of (N/0)r f ratios in the 
artificial composite nebulae from the N/O - O/H relation (the linear 
fit discussed previously) as a function of the difference between the 
oxygen abundance derived by the T e method from global spectra 
and that obtained as the H/3-luminosity-weighted mean abundances 
of components for case B. 

By comparison of Fig. [4] and Fig. |7J it is evident that arti- 
ficial composite nebulae reproduce the anti-correlation between 
(log(N/0)r f - log(N/0)« £L ) and (log(0/H) re - log(0/H) OJV ) ob- 
served in our sample of SDSS galaxies. This can be considered 
in favour of our hypothesis that the enhanced (N/0)r e ratio derived 
for some SDSS objects can be due to the fact that these objects 
are composite nebulae, with a number of H n regions with different 
physical properties contributing to the global spectrum. Hence, the 
electron temperature determined from the auroral-to-nebular lines 
ratio is overestimated, and the oxygen abundance is underestimated 
in such a nebula, i.e., it appears to be a Peimbert temperature fluc- 
tuation effect (see the introduction). 

The middle and lower panels of Fig. [6] show that both the 
ON and NS calibration-based abundances in the artificial compos- 
ite nebulae follow, in general, the abundances of the components. 
However, both case A and B show a gap near 12+log(0/H) ~ 8.1. 
This gap is due to a transition between two regimes in the calibra- 
tions, i.e., by the fact that both "warm" and "hot" Hn regions (ac- 
cording to the classification in Pilvu gin et al. Il20 10b)) make con- 
tributions to the global spectrum of composite nebulae. Therefore, 
such composite H n regions are neither purely warm nor purely hot. 
But distinct c alibration relations for warm and hot H n regions wer e 
in fact used ( |Pilvugineta"ni2010bl : IPilvugin & Mattsson I l2oTTh . 
Hence, there appears to be a problem regarding how these cali- 
bration relations should be applied to such nebulae. It should also 
be noted that the criterion for distinguishing between warm and hot 
H ii regions is somewhat arbitrary. 

We have performed a variety of Monte Carlo simulations, us- 
ing different sets of components and different A(log(0/H)) inter- 
vals, to produce artificial global spectra. We have found that the 
underestimation of the oxygen abundance with the T e method - 
and, consequently, the shift of the position in the N/O - O/H dia- 
gram of the artificial composite nebula with 7Vbased abundances 
relative to its position with mean oxygen and nitrogen abundances 
- is larger for nebulae where the components have large differences 
in physical properties. In particular, the value of the shift increases 
with increasing A(log(0/H). Thus, we have reached the following 
conclusions: 

• If the composite nebula consists of H n regions with similar 
physical properties or a single H n region which makes a dominant 
contribution to the global spectrum, then the oxygen and nitrogen 
abundances derived with the T e method and the ON and NS 
calibrations are in satisfactory agreement with each other, and near 
the mean HB luminosity weighted oxygen and nitrogen abundances 
of the components in the composite nebula. 

• If H ii regions with different physical properties make compa- 
rable contributions to the global spectrum of the composite nebula, 
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Figure 7. Left upper panel. The grey (light-blue) filled circles show the difference between nitrogen-to-oxygen ratios (N/0)r e and (N/0)„, OT „ in artificial 
composite nebulae as a function of the difference between oxygen abundances (O/H)^ - (OfH) mea „ foe case A. The dark (black) filled circles show the 
components. Right upper panel. The grey (light-blue) filled circles show the deviation of (N/0)r e ratios in artificial composite nebulae from the expected N/O 
- O/H relation as a function of the difference between oxygen abundances derived with the T e method from global spectra and those obtained through ON 
calibration. Lower panels show the same but for case B. (A color version of this figure is available in the online version.) 



then the T e - based oxygen abundance can be underestimated. 
Hence, the position of such a nebula in the N/O - O/H diagram 
will be shifted towards lower oxygen abundances, mimicking an 
enhancement of the N/O ratio. 

• The ON and NS calibrations give oxygen and nitrogen abun- 
dances in the composite nebulae which agree with the mean 
luminosity-weighted abundances of their components to within 
-0.2 dex. 

It should be noted that the high-metallicity (12+log(0/H) > 
8.5) calibration Hn regions with measured electron temperatures 
are very few. Thus, they cannot be simulated and we cannot inves- 
tigate how similar (or dissimilar) their abundances estimated using 
the T e method and the strong-line calibrations are. Furthermore, our 
approach does not allow to take into account the contribution of a 
diffuse emission. 



4 DISCUSSION 

Our conclusion that the global abundances derived for composite 
nebulae with the classic T e method may be subject to systematic 



errors is n ot surprisi n g. Fir st, an effect of this kind has been pre- 
dicted by Peimbert (1967) for Hn regions with a n on-uniform 
electron temperature. Second, Kobuln icky et al. I ( 1999) have com- 
pared physical conditions derived from spectroscopy of individual 
H ii regions with those obtained from global spectroscopy for sev- 
eral galaxies with 12 + log(0/H)~ 8.2, and where the auroral line 
[O m]/14363 is detected in both individual and global spectra. They 
found that the standard T e method using global spectra give sys- 
tematic errors in that the oxygen abundances derived from global 
spectra are systematically 0.05 - 0.2 dex below the median values 
computed from spectra obtained in smaller apertures. 

The global oxygen abundances estimated with the 
one-dimensional R2 3 ca libr ation have been discu ss ed b y 
iKobulnickv et al. I jl999h : iMoustakas & Kennicuttl J2006l) . 
while those obtained with the two d imensional R23 calibration (P 
method ) have been discussed by IPilvugin. Contini. & Vf lchez 
d2004j) . It has been found that the abundance inferred from 
the integrated emission-line spectrum of a galaxy using one- or 
two-dimensional ^23 calibrations is representative of the gas-phase 
oxygen abundance at a fixed galactocentric radius (equal to 0.4 
times the isophotal radius), as determined by the abundance 
gradient derived from Hn regions. Here, we have found that the 
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Figure 8. The flux ratio N7 //? 2 as a function of electron temperature f 3 . The 
grey (ligh-blue) filled circles show the calibration Hn regions. The SDSS 
objects with large (> 0.3 dex) deviations from the expected N/O - O/H rela- 
tion are shown by dark (black) plus signs. The Green Pea galaxies with large 
deviations from the expected N/O - O/H relation are shown by open (red) 
circles. (A color version of this figure is available in the online version.) 
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Figure 9. Electron temperature ?3 as a function of oxygen abundance for 
the calibration H 11 regions. 



global abundances derived using the ON and NS calibrations are 
close to the mean abundances. Thus, the calibrations produce 
robust global abundances although both the absolute values of 
abundances in individual Hn regions and of global abundances 
depend on the adopted calibration. 

Our results suggest the 7Vbased abundances in some compos- 
ite nebulae can be incorrect and shift the positions of these nebulae 
in the O/H - N/O diagram towards higher N/O abundances ratio 
or/and towards lower O/H abundances relatively to the positions of 
Hn regions in nearby galaxies. Thus the high (N/0)r e abundances 
ratios obtained in some SDSS and Green Pea galaxies may be in- 
correct because they are composite nebulae. 

Fig. [8] shows the N2/R2 flux ratios as a function of the elec- 
tron temperature for the calibration Hn regions and for the 
SDSS and Green Pea galaxies with large (> 0.3 dex) deviations 
from the expected N/O - O/H relation. In the ON calibration, 
the N2/R2 flux ratios are used as temperature indicators. Fig. [8] 
shows that the measured N2/R2 in SDSS and Green Pea galax- 



ies with f3 ~1.4 have values that are more typical for nebulae 
with ?3 ~1. This means the two electron temperature indicators 
([Oni]/l4959,i5007/[Om]^4363 and N 2 /R 2 ) are in conflict. Thus, 
to reproduce both the measured electron temperature {3 and the 
N2/R2 flux ratios, the composite nebula should contain a compo- 
nent with a temperature f 3 ~1 as well as one with / 3 > 1.4. The tem- 
perature difference between components can be larger than 5000 
K, in some cases. For example, the hottest Green Pea galaxy in our 
sample (with a large N/O deviation), J012910.15 + 145934.6, has 
?3 = 1.75 while its measured N2/R2 flux ratio is typical for a nebula 
with ?3 ~1. 

Fig. [9] shows /3 as a function of the oxygen abundance 
12+log(0/H) for the calibration Hn regions. We see that a Hn re- 
gion with ?3 ~1 has 12+log(0/H) ~8.3. Can a hot component of 
metallicity 12+log(0/H) ~ 8.3 be responsible for the high (3 value 
(~ 1 .4) in the composite nebula? 

We can expect the points in the t 3 - O/H diagram (Fig. [9J 
to be close to the highest attainable electron temperatures for a 
given 12+log(0/H) because of a selection effect: electron tempera- 
tures are easiest to measure in the hottest H 11 regions. In such case, 
an Hn region of metallicity log(0/H) = 8.3 cannot have an elec- 
tron temperature as high as t 3 ~ 1.4. Thus, a hot component with 
12+log(0/H) ~8.3 cannot be responsible for the high electron tem- 
perature measured. 

According to Fig. [9] a Hn region with t$ ~1.4 has 
12+log(0/H) ~8.0. Furthermore, the N2/R2 flux ratio increases 
with decreasing electron temperature. Can a cold component with 
12+log(0/H) -8.0 be responsible for the high N 2 /R 2 ^ 

The typical N/O ratio in Hn regions with log(0/H) ~ 8.0 is 

log(N/0) 1.5. A Hn region with log(0/H) ~ 8.0 and log(N/0) 

= -1.5 would have logA'2/^2 = -0.9 at t 3 ~0.5. However, the line 
fluxes N2 and R2 (normalised to the H/J flux) are one order of mag- 
nitude lower than the measured fluxes. Thus, a cold component of 
metallicity 12+log(0/H) ~8.0 cannot be responsible for the high 
N2/R2 flux ratio measured. 

The measured high N/O ratios and high electron temper- 
atures in some SDSS and Green Pea galaxies can be repro- 
duced by composite nebulae involving components of different 
metallicities. However, the metallicity difference can be as large 
as a factor 2-3, or even larger in extre me cases. It is believed 
that Green Peas are low-mass galaxies JCardamone et al. I [20091 : 
lAmorin. Perez-Montero & Vflchezll2010l : llzotov et al. Il201ll) . The 
metallicity scatte r of H 11 regions in ty pical dwarf galaxies is known 
to be small (e.g. ICroxall et al. 1 120091. This suggests that the mea- 
sured high N/O ratios not in all the Green Pea galaxies can be re- 
produced by composite nebulae and some Green Pea galaxies may 
not be typical dwarfs galaxies. A large metallicity scatter of the 
Hn regions in a dwarf galaxy can arise if this galaxy has under- 
gone an atypical evolution. For example, it can be the result of the 
merging of two dwarfs of different metallicities, triggering a star- 
burst in both components simultaneously. Infall of non-enriched 
gas onto a well-evolved galaxy with a high N/O gas-phase ratio 
would decrease the oxygen and nitrogen abundances, and would 
result in a galaxy with an enhanced N/O ratio. The evolution of a 
galaxy with selective heavy elements loss through enriched galactic 
winds can also result in enhance d N/O ratios (e.g. IPilvugin 1l 1 9931 : 
lYin. Matteucci, & Vladiioboill) . 
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5 CONCLUSIONS 

We have examined the oxygen and nitrogen abundances derived 
from global emission-line spectra of galaxies with 12 + log(0/H) 
ranging from ~7.5 to ~8.5, based on a sample of 281 SDSS galax- 
ies with measured electron temperatures. The oxygen and nitrogen 
abundances in these galaxies were derived with the T e method as 
well as with two recent strong-line calibrations: the ON and NS 
calibrations. 

For r e -based abundances, the positions of some SDSS galax- 
ies in the O/H - N/O diagram are shifted towards higher N/O ratios 
or/and towards lower O/H abundances relative to the positions of 
Hn regions in nearby galaxies. In case of strong-line calibration- 
based abundances, the SDSS galaxies occupy the same area in the 
O/H - N/O diagram as the H n regions in nearby galaxies. 

The global spectra of galaxies have been Monte Carlo simu- 
lated as a mix of spectra of individual components, based on well- 
observed H ii regions in nearby galaxies. Abundance analysis of the 
artificial composite nebulae yields the following conclusions: 

(i) If the composite nebula consists of H n regions with similar 
physical properties or a single Hn region makes a dominant contri- 
bution to the global spectrum, then the oxygen and nitrogen abun- 
dances derived with the T e method and the ON and NS calibrations 
are in satisfactory agreement with each other and near the mean H/3 
luminosity weighted value of oxygen and nitrogen abundances of 
the individual components of the composite nebula. 

(ii) If H ii regions with different physical properties make com- 
parable contributions to the global spectrum, then the T e - based 
oxygen abundances may be underestimated and the position of such 
a nebula in the N/O - O/H diagram will be shifted towards lower 
oxygen abundances, mimicking an enhancement of the N/O ratio in 
the nebula. This effect is similar to the one discussed by Peimbert 
for H ii regions with small scale temperature fluctuations. For com- 
posite nebulae however, this effect appears to be due to temperature 
variations on large spatial scales, caused by varying temperatures 
in different components of the composite nebula. 

(iii) The ON and NS calibrations give oxygen and nitrogen 
abundances in composite nebulae which agree with the mean 
luminosity-weighted abundances of their components to within 
~0.2 dex. ON- and NS-calibration-based abundances for these neb- 
ulae also show a gap in abundance near 12+log(0/H) ~ 8.1. 

(iv) The high (N/O)^, ratios derived in some Green Pea galaxies 
may be caused by the fact that their SDSS spectra are spectra of 
composite nebulae made up of several components with different 
physical properties (such as metallicity). However, for the hottest 
Green Pea galaxies, which appear to be dwarf galaxies, this ex- 
planation does not seem to be plausible. It would work only if the 
HII regions in these galaxies have a dispersion of abundances much 
larger than that typically found in dwarf galaxies. 
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